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ABSTRACT 


The    IEH      Company's  Continuous   Systems      Modeling  Program 
was    used   to   simulate    the      longitudinal   flight   control    system 
of    the   F/A-18    aircraft.       The   model      is   intended    for   use    in 
investigations   cf     aircraft  response   to     flight    conditions 
which  approach      spin    or      stall   and      is   restricted     to    the 
automatic    flaps      up     (AFU)       flight    mode.         Program    outputs 
includ*    stabilator  deflection,      leading     and   trailing   edge 
flap      positions,      and     cress-axis   interconnect      signals. 
Various    stick    forces,    motion   sensor   inputs,    and   air   pressure 
inputs   were   sinulated  to      produce  transient    control  surface 
responses.      Tfcese      computer  generated      responses   exhibited 
characteristics   corresponding   to      predicted   aircraft   control 
surface    movements. 
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I.        INTRODUCTION 

The    Navy    has    experienced  the      loss    of    numerous   aircraft 
during  recent      years    due      to   unintentional      departure    from 
controlled    flight.         The    increased      ccst    and    complexity    of 
modern    aircraft    utilizing    fly-by-wire    flight   control   systems 
have      placed   a      renewed      emphasis      on    understanding     the 
performance    characteristics  of      these   aircraft,      especially 
near    the    limits    of  the   flight    envelope.      The    purpose   cf   this 
thesis  is    to    investigate    a     method    of    non-real   time  computer 
simulation    cf   the  longitudinal   flight      control   system    of   the 
F/A-18  aircraft.        Future   thesis    researchers   at      the   Naval 
Postgraduate    School    will      complete    corresponding   simulations 
cf   the      lateral   and      directional    control      systems.        The 
combination      cf      these      simulations      with      an      existing 
aerodynamic      simulation    program      will      yield   a      complete 
aircraft    stability  and  control  model.    The    primary   purpose   of 
the    model    is      to   investigate   methods  of      designing  control 
augmentation   systems      which  actively      inhibit   or      prevent 
departures    from   controlled   flight.      Othar    uses   cf   this    model 
would   include      the  capabilities      to    test      new   programmable 
memory  configurations,      to  evaluate      new   components  such   as 
ootimal    observers,      to   simulate      degraded    flight   conditions 


(such  as    a    damaged   flight    control   surface)       and    to    recreate 
flight  conditions  during    post    accident    investigations. 

The    method    which      was   chosen    for   accomplishment      of    the 
objectives      of   this      thesis  was      the   Continuous      Systems 
Modeling      Program      (CSMP)  ,        developed        by      the      IBM 
Company   [R€f.    1].        CSMP    is  a     software  package    designed   to 
simulate    dynamic   systems    described      in   tsrms   of    differential 
eguations    and    block    diagrams   normally    encountered   in   systems 
theory.       CSMP    allows    programming    flexibility   through   the    use 
of    thirty-four      pre-programmed  functional    blocks      which   are 
similar   tc    FORTRAN   subroutines.      These   blocks   provide   rapid 
slccqss  to      mathematical    functions,        switching    functions, 
signal  sources,    logic   functions    and    FORTRAN    functions.    Since 
-his   thesis      represented    the   first      attempt  at      the   Naval 
Postgraduate    Schcol    to   accurately      model   the    flight   control 
system  of    a   modern,       highly     augmented    tactical   aircraft   it 
was    deemed    important    to    concentrate      on   the    physical   systems 
rather  than    become   involved   in  the   complexities   of   numerical 
analysis. 

Alternatives   to      CSMP    which      were    considered      included 
analcg      programming      and    FORTRAN      programming.        Analog 
programming      was     not      selected      because      it      is      less 
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representative   of  the   systems   to    be      modeled   and  it    is    less 
accurate    than   CSMP.       The    concept    of    programming    directly    in 
FORTRAN    was      carefully   considered.      Since    FORTRAN      is    the 
source   language      for    CSMP   the      capabilities   of   CSMP      are   a 
subset   of    the    capabilities  of   FORTRAN   itself.       Additionally, 
CSMP    has    restrictions  on    the   number   of   allowable   statements, 
constants,    variables    and    ether  parameters.         Unlike  FORTRAN, 
when    functional    blocks   in   CSMP  are    used,      the    programmer   has 
no    direct    control  of    mathematical   operations   internal    to    the 
functions.      Tie    primary    reasons    for   which    CSMP    was   selected 
were    its    simplified    input    statements,      output   statements   and 
orogram    control   statements   which      facilitated   rapid  program 
writing   and    testing.         Additionally,    the    automatic  time   and 
amplitude    scaling,    data    formatting,      and   compatibility    with 
graphic   display   devices    which   CSMP      provides    are  well    suited 
for    prototype    program  development. 
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II.       METHODOLOGY 

A.        PROGRAM    ORGANIZATION    AND    NOMENCLATURE 

A  detailed  description  of  the  flight  control  laws  of  the 
F/A-18  aircraft  is  given  in  the  McDonnell  Aircraft  Company's 
system      design    report   [Hef-   2].  The   computer      program 

developed    in    this  thesis    is   based      upon   figure    16.1   of    this 
report,    entitled   F/A-18    Longitudinal,    Mechanical,    CAS,       and 
DEL    Control      Law    Mechanization.       This    figure      contains    six 
pages  of    block   diagrams    depicting   generation    of   longitudinal 
control      signals     which      are   valid      for     all      aircraft 
configurations  and      failure  modes.      The  version      of  flight 
control   program    incorporated   in      the   flight   control  computer 
programmable    read     only    memory       (PROM)      utilized      in   this 
simulation    is    8.2.1,     dated   August   31,     1982. 

A  brief   discussion   of   the      F/A-18    flight   control   system 
is   necessary      to    facilitate  a      discussion    of      control    law 
modeling.      All   computations  of  control   laws   are    accomplished 
independently   by     four  channels      of   digital      computation. 
Primary    control    in  the   pitch     axis   is    provided   by  symmetric 
deflection    of   the  horizontal   stabilators.      Full    span   leading 
edge    flaps   and   trailing   edge      flaps    are  scheduled   to   provide 
maximum    lift   to    drag      ratio  during    maneuvering,      high    angle 
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of    attack      and   cruise      configurations.        Roll     control    is 
accomplished   by      ailerons,      differential      srabilators    and 
differential    leading      and   trailing      edge    flap      deflection. 
Directional   control        is    maintained   by    dual      rudders.       The 
thesis  compu-er    program      which  is   contained  in      Appendix    A 
simulates    the   output    of    one     channel   of  digital   computation 
and    calculates      the    angular  oositions   of      the   stabilators, 
leading   edge    flaps  and   trailing    edge    flaps.      This    simulation 
does    not    calculate  rudder     or   aileron   deflection,      however, 
all   electrical   signals  reguired    for   cross    axis    interconnects 
are    provided. 

The    zask    of      programming   the    information    given      in    the 
longitudinal      control      law   mechanization      schematic      was 
simplified    by   two  means.     First,    the    program    restrictions   and 
assumptions    -:c   conditions    of   flight      which    are    discussed   in 
part   B   were   applied.      As    a     result    sections   of    figure    16.1 
which   apply   tc    mechanical   control   laws  and    spin    modes,       for 
example,    were    deleted.      This   reduced   the   number   of   schematic 
blocks  to    be    modeled    by    approximately   one    third.      The    second 
simplification   arose      through  a    system   of      nomenclature    in 
which   nine   control   paths    were   defined     in    order   to   limit    the 
number  of    input      and    output  signals    for     any   specific    path. 
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When    combined,      *hese   nine   paths    fcrm    the    total    longitudinal 
control   law   mechanization. 

Figure    2.1      is  the  overall   longitudinal      control  signal 
block   diagram    which   was    used     in    this   simulation.        It    was 
derived   by  applying    all      program    assumptions   cr    restrictions 
outlined    in      part  B    to   figure      16.1    of    the      system   design 
report. 

Block      diagrams   depicting     the  logic      development   of      the 
component    paths    are    included   in    Appendix   A.    The    nomenclature 
for   each    control    path   which  is      given    in   Table    I   is  peculiar 
to    this    simulation  program. 

Table   II    lists   nomenclature  for      control   signal    groups    which 
are    common      tc    both    the      McDonnell    schematic    and      to   this 
simulation.      Signals      with   common      prefixes    are      numbered 
consscutively    in   the    feed      forward    direction.      The   primary 
feed    forward   input      to   the  CSMP    simulation      is    pilot   stick 
force.      Feedback    signals   include    oitch  rate,    roll   rate,    yaw 
rate,   angle   of   attack,      normal  acceleration  and    differential 
control    surface    commands. 

B.       PROGRAM    RESTRICTIONS     AND    ASSUMPTIONS 

The    task    of   obtaining   a  working   computer    program   to   meet 
the    thesis   objectives   did    not      reguire   the    simulation    of    all 
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Figure  2.1.   Lcugitudinal  Control  Laws 


15 


TABLE    I 
Ccntrcl  Signal   Paths 

PATH  SIGNAL    PREFIX 

Main    path    (open    Icod)  MP 

Pilot   stick    input  pa.h  PS 

Pitch   rate    gyro    path  PR 

Angle   of  attack    sensor    path  AA 

Normal    accelerometer  NZ 

Forward   integrator  path  FI 

Statilatorpath  ST 

leading   edge    flap  path  LE 

Trailing  edge    flap  path  TE 

TABLE    II 
Common   Signal  Nomenclature 

SIGNAL     (OR    PREFIX)  DESCRIPTION 

ALPHAS  Angle    of    attack    (computed) 

ALPHAT  Angle   of    attack    (true) 

F  Function     (air    data   schedule) 

PK  Air  data    schedule    gain    parameter 

PS  Static  pressure 

•    PV  Pitch  axis    storage    location 

QC  Compressible   dynamic    pressure 

RI  Pressure    guotient    (=QC/PS) 

possible    conditions    of  flight.       Thus,      an    assumed   aircraft 
configuration      and      flight      condition        led      to      many 
simplifications    in      the    model.      Each    major      assumption    is 
discussed    below.        Should   future      researchers   desire      to 
construct    a    mere   general      model,      additional    program   logic 
paths  could      be    included    in     this   CSMP     simulation   without 
reguiring    major    revisions    to   the    program. 
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Assumptions: 

1.  The    aircraft    flight  control   electronics   set    (FCES)     is 
operating    normally.       Primary   components  of   this    set   include 
flight   control    computers,      pitch,      roll   and    yaw    rate   gyros, 
accelerometer   assemblies    and    flight    control    panels.    The    FCES 
contains      logic      seguences  for      failure      detection      and 
corrective    control      law    i  nplementation   which      were   ommitted 
from    the    thesis   computer    program, 

2.  The  aircraft  is  under  inner  loop  control.  In  this 
mode  pilot  stick  force  forms  the  primary  control  input. 
Autopilot  functions  such  as  heading  hold,  roll  or  pitch 
attitude    held,    and  altitude  hold   are    not  in   operation. 

3.  The      pitch  control   augmentation   system       (CAS)       is 
providing    longitudinal  control.     Backut)   control   systems    which 
are   net    in    opera-ion    include     the   mechanical   backup  mode   and 
the    direct    electric    link     (EEL)       mode   which    provides  an    open 
loop    signal      frcm  the  pilct      stick    position   sensor      to    the 
stabilator    servoactuator . 

4.  The   aircraft    is   operating   in    the    up   and   away   flight 
phase  in    the   auto  flaps    up    (AFU)     configuration.         This    phase 
reguires    that   the  flap   switch   in      the   cockpit   be   in  the    AUTO 
position   or   that   the    calibrated    airspeed   be   greater  than    243 
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knots.      In    the      AFU    mode    gain   schedules      ar=    optimized    for 
combat   maneuvering   characteristics    in   the    low   to    mid    dynamic 
pressure    regime.         Addit  icnally ,    a    trim  integrator   controls 
load    factor   and,        at   angles   of    attack      above   22    dsgrees, 
proportional   ncse  down  commands    are    introduced. 

5.  The   aain      switch    is  in  the      normal    position.       The 
override    position  of    this   switch      causes    flight    computers   to 
use    fixed    values    for    air      data   schedules  and   a    predetermined 
angle   of    attack      in    control  law    computations.         The  normal 
oosition    of   the      gain  switch   allcws   measured      values   :f   air 
data    ar.d    angle   of  attack    tc  be   used. 

6.  The    aircraft    is  operating    with   weight    off   the  wheels, 
soeed   brake   in,    and    with    cc  external   stores. 

7.  Anti    spin    functions  are  not    simulated. 
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III.       CONTROL    LAW    MODELING 

A.        AIR    DATA    SCHEDULES 

Air   data      schedules    ars      functions   of      sxatic    pressure 
(PS),      compressible      dynanic    pressure      (QC)  ,      and     other 
parameters   such    as  normal   acceleration,      angle   cf    attack    or 
condition    cf    external   stores.      The   longitudinal    control    laws 
contain    inputs    from       13    different   air    data      schedules    which 
are    listed    in    Table    III. 

TABLE   III 
Air    Data    Schedules 

DESCRIPTION     (INPUTS) 

Pitch    forward    loop    integrator    aain   schedule     (RI,PS) 

Fader   on    supersonic  compensation    (QC) 

Srall    margin    on   pitch    forward   loop   integrator    (AOA) 

Trailing    edge    flap   schedule    (AOA.  RI) 

Trailing    edge    flap   schedule    (QC    limit) 

Leading  edge    fiaD    schedule    f&OA.RI) 

Leading   edge    flap   schedule     (QC    limit) 

Leading  edae    flap    schedule     (RI    limit) 

Longitudinal    forward    loop   gain    schedule    (QC) 

NZ    limit    on   AOA   feedback    (NZA) 

Pitch    rate   feedback  gain    schedule     (PS,QC,RI) 

Pitch   rate   feedback  gain   schedule    (QC) 

Longitudinal   inertia!    gain    schedule    (QC) 

Each   schedule      controls   an  output      gain  for      a    particular 
purpose.    Function  32A,   for  example,    yields    a   uniform   initial 
oitch   acceleration  in  response   tc      sharp  inpurs,      with    gain 
decreasing   a-    high  values      of   compressible    dynamic   pressure. 
Other   functions    such    as    function      29U    are    used    to    determine 
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FUNCTION 

F12 

F22 

F23 

F2U 

F25 

F27 

F28 

F29U 

F32A 

F37 

F40 

F68 

F107 

an    upper   cr   loser  limit    to     an  input    signal.      Simulation    of 
these   functions   in  the   CSME      format    required   that    individual 
computer    programs  be    written   to    test  each    air   data    schedule. 
In    all      cases   the      final    result    was      given    in      the    flight 
control   system    design   report,         usually  in    graphical    form. 
Conversion   to   CSMP     format   was   done    by      extraction   of    data 
ooints  from   graphs,       use    of  logic   flow   charts   when   available 
and    direct    employment   of      those    mathematical    formulae    which 
were    given.      Cnce      programmed   in    CSMP,      the      tabular    and 
graphical    output      data  was  compared      to  data    given      in    the 
design   report   for  the  same  test    parameters.    Thus,    the   output 
of    each      function  was      independently   verified      before   the 
function      was   included      in  the      longitudinal   control      law. 
si  mulation . 

The    simulation  is      written  tc   take   advantage      of   CSHP's 
SORT   and    NCSOBT   capabilities   where    acpropriate.         The    NOSCRT 
~oticn   of    CSMP   is  used      for  conditional  logic   and    branching. 
This    statement      allows   the     user    flexibility      in   creating 
sections   of   the   program    in     which  ordinary    FORTRAN   rules   can 
be   used.         In   NOSORT    sections   intermediate      variables    were 
defined,      comparisons      made  and    appropriate      branching    was 
executed.      The      program    statements   were   returned      to    SORT 
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format   as    soon   as  conditional      logic   was    no   longer   required. 
An   early    problem    revealed   that   if      two   or    more    SORT  sections 
are    separated    ty     a    NOSOR1   section   information      will   not    be 
passed  between      the    separated   SORT   sections.         To    prevent 
errors  resulting     from   this  restriction  from      occuring    the 
following      decisions      regarding    the      order      of      program 
statements    were    made. 

1.  The    number  of    NOSORT  sections   should    be    minimized. 

2.  NOSORT      sections      should      be      located      close 
together , allowing  fewer    and  larger    SORT  sections. 

3.  Macros    should  be      utilized    when   possible.      Maoros, 
which  are    similar  to    subroutines,        are  discussed   in   part   D 
under   freguency    averagers. 

The    task    of    arranging    program   statements    to    minimize   the 
number  of    NOSORT    sections   was   greatly    simplified   through    the 
use    of   the      "output    variable   sequence"    tabulation      which    is 
oroduced    as    oart      of    the   CSMP   standard     output.      The   final 
number  of    NOSORT   statements  in   the    computer   simulation    could 
have    been    further  reduced    by      this    method,       however,      this 
would   have   reguired   the    movement    of    large    program   blocks    and 
caused  a      deviation    from    the     logic   oath   used     in   program 
development.      For  example,   the   NOSORT    sections   of    functions 
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12    and      40    could    have      been   combined.         This      would  have 
increased    difficulty    in    program    debugging      and   mads   areas   of 
the    program    which  require     combinational   logic   less  apparent 
to    readers. 

Appendix   C      contains    documentation     for   the      air   data 
schedules.      Each   schedule    is   included   in   a   complete  computer 
orcgram    which    produces  tabular  and      graphical   data   to    match 
the    characteristics    of      figures    in   the   design      report.       In 
incorporating   these    air    data    schedules  into    the    longitudinal 
control   law   simulation  all  schedules      exceot    for   function    24 
were    placed      before    the      main    computational      body   of      the 
program.       Function  24,      the     trailing   edge    flap   schedule, 
reguires    conditional      logic  and    contains    computed      angle   of 
attack  as      an    input.      Sines      computed   angle   of      attack    is 
generated    in   the    normal    accelero  merer    path,      function    24    was 
placed     immediately    before     the    trailing      edge    flap     path 
computations . 

B.        FADERS,       AKAICG    TO      DIGITAL    CONVERTERS,       MECHANICAL 

EREAKCUTS 

The   purocse  of   faders  is   to  eliminate   large 
discontinuities  in  gain,  permitting  gradual  change  in  output 
from  old  values  to  new  new  values  at  a  desired  sample  rate. 
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Faders  are    located   in    Figure    16.1    at   the    outputs   of   air    data 
schedules      and   gain       schedules   which      are    dependent      upon 
aircraft    configuration.         Discontinuities    may  arise     as    a 
result  of    a    change   in      physical    measurements    such    as    dynamic 
pressure,      change  in    aircraft      configuration    such    as   speed 
brake  extension,     or    change  in      mcde  of  flight.      The   lower 
limit   for    signal    MP4,      for  example,      changes   from    a  gain   of 
-50.0  to    0.0    as    a  result      cf   spin  entry.         Since   the   thesis 
computer      program  assumes      that      electrical   signals     vary 
smoothly    and    is   restricted  to   up    and   away   controlled    flight, 
faders  are    modeled  in   freguency    only. 

Analog   to    digital  conversion      and    frequency    matching   are 
obtained    ty    a   sample    and      held  process.      Sampling   times   are 
generated    by   the    CSMP  functional    block    IMPULS   which  produces 
a   time   series      cf  unit    impulse   functions      with    a  specified 
start  time    and   period.      Since   these    pulse    trains  are  used   in 
several   areas   cf   the    program,      impulse   functions   of  23,      40 
and    80  hertz    are    included    immediately   following   the   air    data 
schedules.        The    zero   order  hold    function    ZHOLD   keeps    signal 
gain   constant    throughout    the   pulse   period. 

d   mechanical    breakout    force    cf      two   pounds   is   modeled   in 
the    pilot    stick   input  path.  When    large   stick   forces    are 
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applied,    the   resultant   signal   is    calculated   with    a  reduction 
in    magnitude    cf   -wo    pounds.   This    modeling    is  accomplished    in 
CSMP    by   the   deadspace    (DEADSP)    functional    block.    Figure    16.  1 
of    the   design    report      depicts   electrical   signal   "deadbands" 
which   are      conceptually    identical   to      mechanical   breakouts. 
Appendix   D      cortains    a      program    which      displays    deadspace 
outputs   in   tatular  and   graphical    form. 

C.        ALIASING    FILTERS     AND    SIGNAL    LIMITERS 

Ihe    longitudinal    control    laws  cf   the    F/A-18    include    five 
aliasing    filters    which  are   modeled   as   first   and   second    order 
Laplace    transforms.         The    first    order     transform   is  of    the 
tvpe    A/(Bs+1)    and  is    in   the   pilot   stick  input   path.    This    lag 
type    filter    with   one    real    pole      is    converted   tc    CSMP   format 
by    the  functional  block    REALPL.     A  reguired    initial   condition 
is    the   value   cf   the    output     signal   when  time    is    zero.      This 
may    be   determined  arbitrarily    by    the      user,      but    was   set    to 
zero   for   this   simulation. 

Second   order      filters   are   present      in    the      pilot   stick 
input  path,      pitch  rate    gyro    oath,      angle    of   attack  sensor 
oath    and    normal      accelerometer   path.      Each    is      of   the    form 
A/  (Bs2  +Cs  +  D)     and   represents  an  underdamoed   system.      The   CSMP 
functional    block      for   complex      poles    (CMPXPL)         is   used. 
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Initial   conditions   are      tke  value    cf   the      output   signal    and 
time    rate    of      change    of    the  output    signal      when   time   equals 
zero.      Natural    frequencies     varied    from   4.34    hertz      in    the 
pilot   stick    input  path  to    33.3      hertz    in    the    angle    of    attack 
sensor  path.      The  appendix   contains   a   computer    program    which 
demonstrates   the   CSMP   outputs   of    both   first  and    second    order 
aliasing    filters.      Unit    step   inputs    wera    introduced  to    each 
filter  to      generate    transient   responses.      The      first   order 
filter  produced   an  exponential  rise    to   steady  state  with   the 
correct    time    constant.      Characteristics  of    the    second    order 
filters   such   as   rise    time,   peak    time,    maximum   overshoot,    and 
settling      time      compared        favorably      with      theoretical 
re  suits   [Eef .    3]. 

Sianal    limiters    restrict   the      maximum    or    minimum   values 
of    an  output   signal.        Stafcilator  surface    deflection,      for 
example,       is   limited    to    10.5   degrees  trailing   edge   down    and 
24    degrees   trailing    edge    up.    The    CSMP    functional   block    LIMIT 
allows      direct   specification      of    lower     and    upper      signal 
limits. 

D.        FREQUENCY    AVERAGERS    AND    RATE    LIMITERS 

A   characteristic    cf   the  F/A-18      flight   control   system   is 
that    various   signal    paths   operate   at    frequencies  of  20,      40 
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or    80   hertz.       When  signals  are   combined   mathematically ,    the 
inputs  are    first    converted   to     a    common  frequency.      Normal 
accelerometer    path      signals  are    computed      at    40      hertz    and 
combined    with   outputs   of    the   forward   integrator   path,       which 
operates    at    20   hertz.    In    this  case   a    20   to   40   hertz  averager 
is    present      between    the    integrator      path   and      the   summing 
junction    connected  to  the      normal  accelerometer   path.      The 
algorithm    used   for  the    20   to    40    hertz    averager    is   based   on   a 
Drocedure    given      in    the      Flight    Control      Electronic  System 
Report  [Ref.    4],        The      averager    was   required      to   generate 
signal  values    at     twice   the  rate    of   the      incoming   pulses    by 
linear  interpolation    between  amplitudes     of    the   two   previous 
signals   at      20   hertz.  The    formulas      used  to      generate 

intermediate   signals    were: 

Z40    =    Z40Z1    ♦    DEL  (3.  1) 


DEI    =    (Z20    -Z2CZ1)    /    2.0  (3.2) 

where  : 


Z40        =   current   value   of    the    40   hertz   signal 
Z40Z1    =   previous    value   of    the   40   hertz   signal 
Z20        =   current  value   of    the   20   hertz    signal 
Z20Z1    =   previous    value   of    the   20  hertz   signal. 
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Conditional    logic   was      used    to    keep   the      output   signal 
equal   to    the    input   signal    at      times    when    ths    20    hertz    impuse 
function    was    egual  to   one.      The    algorithm    was   initialized    by 
letting    "previous  values"    be   equal   to   input    values    when   time 
equals  zero. 

The    presence     of    numerous      frequency    averagers      in    the 
control    laws      would    have    reguired      the    repetition     of    many 
statement    blocks    without      the   use  cf   CSMP      program   MACROS, 
which      are    similar      to      FCETRAN    subroutines.         Frequency 
averagers    in    tie   thesis   computer    simulation    were   included   in 
MACRCS  and    placed  at    the    beginning   of   the    program.      A    MACRO 
may    be  used   several    times    within   a    program.    Input   and   output 
variables    are   given    dummy    names,         yet   the    MACRO  is  invoked 
with      a    unique      name    which     is   assigned     in   a      function 
definition   statement.      A    limitation   to   the    use   of    MACROS    is 
that      variables    which      are     defined      in   MACRO      structure 
statements   are      not    available   for      output    unless      they    are 
designated      as      arguments   in     the      function      definition 
statement.         Additionally,    certain    functional    blocks   such   as 
REALPL,       CMPXPL    and    INTGRL   cannot   be    used    as   arguments    or   as 
parts  of    a    MACEO   structure  statement. 
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The    PROCEDURE  function  of   CS MP    was   used    in    each  MACRO    to 
cause   statements     to    be      executed   in      tha    order      of  their 
appearance.      Each  structure     statement   of    a   MACRO      will    be 
individually    scrted      unless  PROCEDDRE   is      specified.       All 
statements   contained    within  a   PROCEDURE  function   are   treated 
as    a    single   block  which    can      be    moved    but    not   rearranged   by 
the    CSMP    translator. 

FORTRAN   subprograms      were   not   used   in      this    simulation 
because   the   size   cf    the      longitudinal    portion   of   the   program 
did    not   approach   the    maximum   limits   of   CSMP.         Inclusion   of 
the    lateral    and    directional  control      systems   in   this   program 
will    reguire   that      measures  be  taken    to     remain    within    the 
allowable    number   of    structure     statements,      NOSORT   sections 
and    MACROS.      The     CSMP  translator   is   not      used   to   process 
FORTRAN      subprogram    statenents.         Since      the    number     of 
subprogram   statements  is    net   counted  the   overall   size   of   the 
orcgram    may    te   increased.      The  capabilities   of    the   computer 
system  library   may      be   utilized    by    subprograms      through   the 
use    cf  the   CALL    statement.    The   CALL    statement   mustr    however, 
be    included    in    a    NOSORT    or      PROCEDDRE    section.       A    method    to 
invoke  subprograms  exists    which    does      not    reguire   NOSORT   or 
PROCEDURE    sections,       but    it  is  valid      only   for    two   or    more 
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output   variables      and   it    requires      a   specific      format    for 
arguments.         A      final   restriction   to      the    use      cf   FORTRAN 
subprograms   is   that    certain  CSHP      functional   blocks  such    as 
ZHOLD,   IMPULS,    and  CMFXPL    are   not   allowable. 

The   longitudinal     contrcl   laws     contain    rate      limiters 
which  operate    at      frequencies   of    20   hertz    and      80    hertz   to 
restrict    the      speed    of   leading      edge    flap      movement.       The 
algorithm    for      these    limiters  compares     the    value     of    each 
incoming   signal      with  the      value    cf      the    previous     output 
signal.      The      magnitude    of     the    difference      between  these 
signals    is    processed      by    the   LIMIT   functional      block   which 
generates    the   current      output    signal.      The   20      hertz   rate 
limiter,      for    example,      allcws   a      maximum    change  in  output 
signal   value      cf    0.9    degrees     during   each    period      of   0.05 
seconds.     A    listing  cf   the    program   which  was    used  to  test   the 
rate   limiter    is    contained    in     appendix  D.         Correct  outputs 
were    observed      fcr  bo-h      increasing   and      decreasing  input 
signals. 

E.        DIGITAL       FILTERS,       DIGITAL      TO    ANALOG       CONVERTERS, 

SERVOMECHANISMS 

Three      types   of      digital   filters      are      used   in      the 
longitudinal    ccn^rol    laws.      Z   filter   number   P2    is    a  lead-lag 


29 


type   controller   in  the   pitch     rate    gyro   path   which   operates 
at    20     hertz.      Z     filter    number      P9 ,      the      forward   loop 
integrator,      operates   at    20   hertz  and   compares   the   aircraft 
response    to   the   maneuver    command.      The   output   signal   drives 
the    stabilator   servoactuatcr  to    reduce  the   maneuver  error   to 
zero.      This    allows   the   aircraft    to    to    be    automatically    kept 
in    a     hands   off      condition  since      the    forward      integrator 
eliminates    unccm landed  nornal  acceleration.      Z   filter    number 
P8    is  a    structural  notch    filter      which   operates    at   80    hertz. 
It    attenuates    aeroelastic   bending   which      is   detected  by   the 
moticn  sensors. 

In  the    simulation   each   filter      was   developed  in  its   most 
aeneral    form      for  inclusion  in      a    MACRO.  The    following 

equation    for      lead-lag  filter   number      P2    is   given      in    the 
Schematic    Cesign   Report    [ Bef .    5]. 

PR4  (1+PK11*(1-EK12) )Z    -    (PK1 1+1)* (1-PK12) 

PR3  Z     -    (1-PK12)  (3.3) 

It  is  modeled  in  the  thesis  simulation  as 

FOUT  AZ  -  B 


FIN  Z  -  C  (3.4) 

where  A,   B  and  C  are   constants.   The  right  shifting  and 

linearity  properties  of   the  z  transform  are   used  to  solve 

explicitly  for  the  variable  FOOT  fRef.  6]. 
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FODT  A    -    B  (Z-  1) 

FIN  1 .0   -    C (Z-1) 


(3.5) 


Cross   multiplication    and    rearrangement   yields: 


FO0T=    A*FIN    -    B*FIN(Z-1)     ♦    C*FOUT(Z-1) 


(3.6) 


which   is    described   in  the    simulation   as: 


FOUT=    A*FIN    -     B*FINZ1     +    C*FOUTZ1 


(3.7) 


This      method,      which      is      termed    direct      realization 
programming,      was  also   used  in  the    development   of    the    notch 
filter  and      the    forward    loop   integrator.         The    eguivalent 
Laplace   -ransfcrm      for  each      longitudinal    flight      control 
filter  is    listed    in    the      design    report.        This    permitted   a 
cress   check    of      z   filter    performance    which      is   included    in 
appendix    D. 

A   specific   method  of    integration   may   be    specified    in    the 
terminal    pcrtion   of    a   CSME  program.         In   the    case   of  flight 
control   simulation  the  Runge-Kutta      Fixed    Step   Size    (RKSFX) 
method   was   utilized    to  ensure      that    integrations   would    only 
occur  at    the    desired    sampling      rate.      The    highest    sampling 
freguency    in    any    axis  of    the  F/A-18   flight   control   system    is 
80    hertz,         thus   the    CSMP      integration   interval      DELT    was 
soecified    as    0.0125    seconds. 
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The  variable  KEEP  is  used  in  CSMP  to  indicate  that  the 
end  cf  a  valid  integration  step  has  bean  reached.  KEEP  is 
set  equal  to  one  when  this  condition  is  met.  During  trial  or 
intermediate  integration  steps  KEEP  will  equal  zero.  Each 
MACRC  contains  conditional  logic  which  allows  calculations 
to    be   performed    only    when    KEEP  equals   one. 

Conversion   cf  signals    from  analog   to   digital   form    in   the 
F/A-18   occurs   as   the    signal  reaches   the  servo mechanism.      The 
quantizer    functional    block    (QNTZ  B)     is   employed    to    accomplish 
the    analog   to    diqital     simulation.      The   transfer   functions 
used    by      the    stabilator    and      flap   servcmechanisms     are    not 
oublished    ty    the   manufacturer   of    the  aircraft.      In   order    for 
the    thesis      computer    program     to    generate      control   surface 
positions    the    response   characteristics      cf    a   Parker  Hannifin 
fly-by-wire    actuation  system   were   incorporated.    The  selected 
actuators    were      designed    for      use   wi-h      all-digital   flight 
controls    and    are      modeled    as    second   order      systems.         The 
stabilator    transfer    function   is 

X  1600 


Ei  S2*    56S    ♦     1600  (3.8) 

where  X  is  the  actuator  position  in  degrees  and  Ei  is  the 
position  ccmmand.  The  transfer  function  for  both  leading 
and    trailing   edge  flaps    is 
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X  400 

Ei  S*  +    28S    +4  00 


(3.9) 


The    advantages   of  using    the  second      order    model    instead    of   a 
first  order   model  are  that  the  faster   rise    time   more   closely 
represents   the      physical   actuators     and   that      the   natural 
frequency    and   damping   ratio  may    be   independently   modified. 

F.        PROGRAM    TESTING     METHODS    AND    RESULTS 

The   thesis   computer    program    was      tested   on   three   levels. 
The    lowest      level  involved     evaluation   of      the    individual 
signal  blocks   cf     figure    16.1    in    the      system    design  report. 
Sections    A   through  E    of      this   chapter   describe   signal    block 
modeling    techniques      and    appendicies      C  through      G    contain 
testing   programs.        To  obtain   verification    of    proper   program 
operation      it    was      desired  to      create      input   signals      of 
realistic    value   which   would   produce      a   time    varying  output. 
The    rate    limiters,      for    example,      were   tested   with  an    input 
signal  which    rose  exponentially    to      a   limiting   value,      then 
decreased    ex potentially.  Output   signals   at      the   desired 

freguency  were  observed  for  incoming  signals  of  positive  or 
negative  slope  and  for  incoming  signals  which  were  within  or 
beyond   the    rate    limit. 
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The    second    level      of    program    testing    involved      the    nine 
signal   paths    which  are   listed    in      Table    1.       In    most,  cases    a 
s~ep    function    was  used   as      the  input   signal.      Intermediate 
signals      were        observed      to        determine      continuity, 
f  reguency  ^ime   constants,        and    conformance      with    limiting 
values. 

The      highest    level      of  program      testing   required      the 
aer.eraticn   of      control  surface      deflections   for      specific 
combinations   of    pilot   stick  force  and   motion   sensor   feedback 
signals.       Neither     aircraft   flight   data   nor      the   McDonnell 
flight  control      simulation   data    were   available      for   direct 
comparison    with    the      outputs   of    the    thesis      program.       The 
description   of   control   characteristics      given   in  the   design 
reoort   was   used      to    make    a  qualitative      analysis   of  program 
operation . 

One    specific   condition  of   flight    was   selected   as   a    basis 
for    comparison   of  control      responses   to  various   combinations 
of    input      signals.      This    "base      condition"    of      flight    was 
defined   such      that  inputs      from    motion     sensors   could      be 
superimposed      on    base      condition      inputs      to    enable     an 
investigation    cf   the    effects   of      each    parameter.        The    base 
condition    for      program   testing  was      selected   to      model    an 
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Figure   3.1.      Base   Condition   Stabilator   Response 


aircraft   operating  at   20,CC0   feet   and    250    knots.      This    fixed 

the      values   of      static  and     dynamic      pressures    for      each 

simulation    and      permitted    nanual    verification   of      the    gains 

produced    by   the    air    data    schedules.         Additionally,      a    step 

function    representing      six   pounds      cf    force      in   the     aft 

(positive)       direction  on    the   control      stick   was    applied   at 

time    0.0    seconds.   All   moticn   sensor   inputs    were    held  at    zero 

so    that   their    effects   could      be    individually    studied.       The 

initial   deflections    of  the     stabilator,      leading  edge    flaps 

and   trailing    edge      flaps    were   set   tc   zero      degrees   and    all 

initial   conditions   for   filters   were   set   to    zero. 
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Figure    3.2.      Base    Condition  Leading  Edge   Flap   Response 

Figure   3.1    depicts     the   movement   of   the      stabilator    in 
response      to   tke     base   flight      conditions.       A      transient 
oscillation    is   produced      in  the    first   second      which  results 
orimarily    from   the      second   order    filter  in      the    pilot    stick 
dynamics    path.         After  this     oscillation   has      decayed   the 
stabilator    continues    to    deflect    at      a   nearly   constant    rate, 
since   feedback   inputs   are    supressed.      The    stabilator  reaches 
the    limiting    value  of   24    degrees      trailing   edge   up   after    9.6 
seconds. 
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Figure    3.3.      Leading  Edge  Flap    Response    to    AOA   Feedback 


Figure    3.2   is     a    plot    of   leading      edge    flap    deflection 

(LEFLA?)     versus   time    for    the    base    flight    condition.    A    steady 

state   flap    deflection  of    4.8     degrees   is   achieved   after    1.5 

seconds.     Since    LEFLAP  is    a   function    only    of   angle    of   attack, 

static   pressure      and    dynamic   pressure      it    was      desired    to 

observe   the   variation  in    LEFLAP    with   angle    of   attack.     A   ramp 

type   increase   in     angle    of  attack   sensor      input    (AA1)       was 

superimposed    upon      the  base   flight   condition      beginning   at 

time    1.4    seconds. The    resulting  schedule  of    leading   edge    flap 

deflection    is    shewn    in  Figure    3.3   and    is   consistent   with   the 

functional    description   givan   in    the   design    report. 
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Figure    3.4.      Base    Condition  Trailing  Edge   Flap  Response 

Trailing    edge      flap    deflection      was    modeled      similarly    in 
Figure   3.U    for   the  base    flight      condition    and   in   Figure    3.5 
for   the    condition  in    which   AA1      is    a   ramp    function   starting 
at    time    1.4    seconds. 

lo  determine    the    effect  of  forward    pressure    on   the    pilot 
siick   a   step   input   of  -12. C   pounds   was   superimposed  upon   the 
base    flight   condition   at    time    4.0   seconds,       which    simulated 
an    instant anecus   reversal   cf   stick   force.      Figure    3.6    shews 
that    the    direction  of  statilator   deflection   changed  abruptly 
upon    introduction  of      the    new    stick    force,         and  that    the 
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Figure    3.5.      Trailing   Edge   Flap   Response    to    AOA 'Feedback 


magnitude    cf   the   steady    state      rate    of   stabilator    deflection 
was    approximately  unchanged.      Pilot      stick    force   inputs   of 
larger  magnitude    were   simulated,       but      the    results   are    not 
shown.      In    these    cases   the      rates   cf    stabilator    deflection 
increased      while    initial      oscillatory    behavior      exhibited 
characteristics   similar   tc  base    condition    response. 

The    influence  of    angle  of      attack    feedback    on   stabilator 
position    is   displayed   in    Figure      3.7.      As    described  in    the 
design    report,    angles  of    attack    in   excess    of    22    degrees    will 
generate    a    proportional    stabilator      command   causing  the    nose 
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Figure    3.6.      Stabilator   Response   to   Stick   Force   Reversal 

of    the  aircraft   to  pitch    down.       In   this   simulation   the    angle 
cf    attack    input      function    AA1    was   set   equal      to   22   degrees 
plus   a  ramp   ~ype    increase   of     one   degree    oer    second  starting 
a-    time    4.0    seconds.    The    resulting   stabilator    deflection   was 
identical    to      -hat    for   the     base    condition    until      time    4.0 
seconds   due   to   the  action    of  the      -22    degree    bias    in   the    AOA 
feedback    path.    The  nose    dcwn   pitching   effect    of   AOA   feedback 
is    apparent    at      all    later   times.      In    response      to    the    base 
condition    the    stabilator    had   reached    its   maximum   limit    of    24 
degrees    trailing      edge   down      at    time      9.5    seconds.        The 
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application   of   AOA   feedback  restricted   stabila-cr   deflection 
to    20.9    degrees   at  time    9.5  seconds. 


Figure  3.7.       Ingle  of   Attack   Damping 

Normal  acceleration  damping  is  shown  in  Figure  3.8.  A 
ramp  increase  in  normal  acceleration  equal  to  1.0  g's  per 
second  was  superimposed  on  the  base  condition  at  time  4.0 
seconds.  Ths  output  signal  varies  smoothly  due  to  the  fact 
that  normal  acceleration  path  outputs  are  processed  by  the 
forward  integrator  path.  The  stabilator  reverses  its 
direction    of    movement  within   one    second  of    the    time  feedback 
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is  introduced.  In  this  cas<=  the  stabila-or  reached  its  limit 
of  deflection  of  10.5  degrees  trailing  sdge  down  a+  time  9.0 
seconds. 


TIME    (seconds) 
Figure    3.8.      Normal  Acceleration   Damping 

According      to  the      design      reocrt,      the      predominant 
contribution    to    pitch   damping   is    generated    by   the   pitch   rate 
gyro    path.      This     is    because   the    pitch   rate      signal  PR5   is 
summed   with      the   main   path     signal    MP9   downstream     of   the 
forward    integrator  as   shown   in   Figure   I.       Figure   3.9   depicts 
stabilator   position    versus  time    for   a    ramp   increase  in    pitch 
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rate,  beginning  at  time  4.0  seconds.  The  change  in  direction 
of  stabilator  icovemen-  is  much  more  rapid  than  that  produced 
by    angle    of   attack   or  normal  acceleration   damping. 


0.00  1.00  2.00  3.00  11.00  5.00  3.00  7.00  8.00  9.00  10.00 


TIME    (seconds) 
Figure    3.9.      Pitch    Rate   Damping 

The   thesis   computer    prcgram      used   a   significant  portion 
of    allowable   CSMP  program    size.    Table    IV,    which    is   extracted 
from    the    "translation   table"    section   of  CSMP   output,      lists 
the    areas   of   the    program      which    most   closely    approached  the 
size    limits   of  CSMP. 
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TABLE    IV 
Program    Size    Restrictions 


Current 

Parameter  Program 

MACRO    and    statement    outputs  206 

Statement    input    work    area  422 

Parameters-function   generators  43 

History   and    memory   block   names  21 

MACRO    statement    storace  85 

SORT    sections  5 

Maximum   statements   in   section  171 


CSMP 

Maximum 

600 

1900 

400 

50 

125 

20 

600 

IV.       CONCLUSIONS     AND    RECOMMENDATIONS 

Ihe    computer    program    developed      in   this    thesis    simulates 
the    performance    cf  the   longitudinal   flight   control   system   of 
the      F/A-18      aircraft      by     gensrating     control      surface 
deflections   and    cross  axis  electrical   signals.    The   responses 
to      flOA,      normal     acceleration   and      pitch    rate      feedback 
correspond      to   the      descriptions    given      by    the      aircraft 
manufacturer . 

The    use   of   CSMP    simplified      the    task    of   system    modeling 
by    provision      cf    pre-programmed      functional  blocks      and   a 
flexible      format    for      outputs.      The      CSMP    orogram     size 
restrictions    dc      not    appear     to    be      a    factor      which  would 
prohibit    the   addition  of    the     lateral   and    directional    flight 
control    systems   to   this    simulation.      It   is    recommended    tha-1: 
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future  expansion    cf   this    program    be   dona    using   techniques   to 
conserve    program      size.       FORTRAN    subroutines      and    computer 
library   functions      should    be   utilized      dua    to      the   limited 
MACRO  and    NOSORT   capabilities   of    the   CSMP    translator. 
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APPENDIX    A 


FLIGHT   CONTROL    SYSTEM    COMPONENT    PATHS 


PILOT    STICK    INPUT    PATH 


/T>ILOT 
/    FORCE 

STloT\ 
INPUT     > 

| 

FAILURE 
DETECTOR    (NOT 

1             (Lba; 

1 

MODELED) 

J 

PS1 

I 

PS5 

MECHANICAL 

BREAKOUT    (2  LB 

NIN) 

0EA0BANO    (-/- 
0.06    IN) 

} 

PS2 

r 

PS6 

STICK    DYNAMICS 
(SECOND   OfiOER) 

STICK   GRAOIENT 
(FUNCTION  20) 

\ 

PS3 

I 

i 

PS7 

r 

ALIASING  FILTER 
(FIRST   ORDER) 

LIMITER 

j 

PS4 

! 

] 

PSB 

1 

ANALOG  TO 

OICITAL 

CONVERTER   (BO 

HZ) 

PS5 

DIGITAL  FILTER 
P5 

' 

1 

V 


U6 


c 


PITCH  RATE 
CYRO 


PR1 


ALIASING  FILTER 
(SECOND  ORDER) 


PR? 


ANALOG  TO 

OICITAL 
CONVERTER 


PR3 


FAILURE 

DETECTOR   (NOT 

MODELED) 


PR) 


OICITAL  FILTER 

Pi 

(N/A) 


PITCH    RATE    GYRO   PATH 


PR4 


LIMITER 


RV4 


ia 


6 


20  HERTZ 
SANPLE/HOLO 


I PV4A  I 


(F6S   »  F32A)  ♦ 
F40 
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ANGLE  OF  ATTACK  SENSOR  PATH 


ANGLE  OF 
ATTACK 
SENSOR 


AA5   (ALPHAS) 


AA1 


ALIASING  FILTER 
(SECOND  ORDER) 


AA2 


ANALOG  TO 

DIGITAL 
CONVERTER 


AA3 


FAILURE 

•DETECTOR   (NOT 

NOOELED) 


AA3 


ANGLE   OF   ATTACK 
COMPUTATIONS 


< 


TO  LEADING 

AND 

TRAILING 

EDGE  FLAP 

P»THS 


-22.0  OEGREE 
HAS 


NZA    (FRON 
NORMAL 

C*}>A?H?TEP 


ALPNAT   (TO 

LATERAL  AND 

DIRECTIONAL 

SYSTEMS) 


PV1 


ro  MAIN  PATH" 


?0   HERTZ 
SAW>LE/HOLD 


PV1A 


TO  FORWARD 
INTEGRATOR 
PATH 


U8 


NORMAL  ACCELEROMETER  PATH 


f      NORMAL       \ 
(  ACCELEROMETER 


V 


AA1 


ALIASING  FILTER 
(SECONO  ORDER) 


AA2 


ANALOG  TO 

DIGITAL 
CONVERTER 


AA3 


FAILURE 

DETECTOR    (NOT 

MODELED ) 


AA3   (NZA) 


-& 


N24 


"RRI    *  RR1   • 
6.l5?9E-« 

(LATERAL 
STSTEN) 


TO  LATERAL 
SYSTEM. 


DIGITAL   FILTER 
P5 


NZ5 


LIMITER 


TO  ANCLE   OF 

ATTACK 
SENSOR  PATH 


PV2 


?0  HERTZ 
SAMPLE /HOLO 


L3 


F32A  *  3.5 


U9 


FORWARD    INTEGRATOR 
PATH 


POSITIVE   LIHITEB 


FI4 


FI5 


DIGITAL   FILTER 
P9 


VARIABLE   LIMITER 


FI« 


f ALPHAS 
/   (FROM  AOA 
I  SENSOR  PATH) 
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STABILATOR  PATH 


DSl 


INFERENTIAL 

STABILATOR 

COMMAND 


OS1 


ST6 


LIMITER 


ST3 


LIMITER 


ST7 


DIGITAL   TO 

ANALOG 
CONVERTYER 


OIGITAL    TO 
ANALOG  CONVERTER 


ST4 


STt 


ST5 


LEFT 
STABILATOR 

SERVO- 
MECHANISM 


RIGHT 
STABILATOR 

SERVO- 
MECHANISM 
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LEADING  EDGE 
FLAP  PATH 


5   (FROM  ^\ 
1  SENSOR      > 
PATH) 

1 

F29 

1 

' 

] 

' 

DIGITAL   FILTER 
Pll 

VARIABLE  LIMITER 

' 

LE2 

LE3 

20  HERTZ 
SAMPLE/HOLD 

RATE  LIMITER 
(20  HERTZ) 

LE2T 

LE4 

,     ' 

F27 

> 

A 

1' 

\T 

20   HERTZ  TO  SO 
HERTZ  AVERAGER 

LI2T1 

» 

VARIABLE  LIMITER 

F28 

LE2T2 

1 

' 

1 

1 

K 

•> 
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LEADING  EDGE    FLAP 
PATH    (CONT'D) 


leu 


limited 


LEI? 


18  DEGREE/SEC 

RATE   LlKlTER 

(80  HERTZ) 


LE13 


LIMITER 


INFERENTIAL^ 

PL AP( LATERAL 
SYSTEM) 


OLE1 


VARIABLE 
DEAOBANO 


DLE2 


0LE3 


LE6 


LIMITER 


LE7 


18  DEGREE/SEC 

RATE  LIMITER 

(80  HERT2) 


LE8 


LIMITER 


OIGITAL  TO 
ANALOG  CONVERTER 


LEU 


OIGITAL    TO 
lANALOC  CONVERTER 


LE9 


LE15 


LE10 


LEFT  L.E. 
FLAP  SERVO- 
MECHANISM 


RIGHT   L.E. 
FLAP  SERVO- 
MECHANISM 
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TRAILING  EDGE   FLAP   PATH 


AA5    (FROM 

AOA   SENSOC 
PATH) 


OICITAL   FILTER 
P12 


TE2 


20  HERTZ 
SAMPLE /HOLD 


TE2T 


IS! 


VARIABLE   LIN ITER 


TE2T1 


6 


5U 


TRAILING  EDGE   FLAP   PATH 
(CONT'D) 


& 


TE7 


LINITER 


DIGITAL   TO 
(ANALOG  CONVERTER 


OTE1 


INFERENTIAL^ 

T.E.   FLAP 

(LATERAL 

STSTEN) 


TE4 


LIMITER 


TEI 


DIGITAL   TO 
ANALOG  CONVERTER 


TE5 


TE» 


TE« 


LEFT  T.E. 

FLAP  SCRVO- 

(CCHANISM 


RIGHT   T.E. 

FLAP  SERVO- 

»CCMANISM 
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ilPENDIX    E 


FLIGHT    CONTROL    COMPUTER    PROGRAM 

4c  **************  ***************  *************************** 

*  **  ** 

*  **  LONGITUDINAL    FLIGHT    CONTROL    LAW    SIMULATION  ** 

*  **  ** 

*  ******************************************************** 
* 

*  *************************************************** 

*  *****INTEGRATING, LAG, AND    NOTCH    DIGITAL    FILTERS***** 

*  *************************************************** 

MACRC    FOUT  =  ZINT(FIN,  KA,  KB  ,  IMP,  FOUTZ  1) 

PROCEDURAL 

IF(IMP.NE.     1.0)     GO    TO    10 

IF(KEEP.NE.  1.0)     GO    TO     10 

IF  (TIME.  EC. 0.0)     GO    TO    10 

FOUT=KA*FIN    +    KB*F0UTZ1 

FCUTZ1=F0UT 
10       CONTINUE 
ENDMAC 
MACRO    FOUT=ZLAG(FIN, KA,KB,KC,IMP ,FINZ1 ,F0UTZ1) 

PROCEDURAI 

IF(IMP.NE.     1.0)     GO    TO     10 

IF  (KEEP.  NE.  1.0)     GC    TO    10 

IF(TIME.EQ.O.O)     GO    TO     10 

FOUT=KA*FIN    -    KB*FINZ1    +    KC*FOUTZ1 

FOUTZ 1=FOUT 

FINZ1=FIN 
10       CONTINUE 
ENDMAC 
MACRC    FOUT=ZNOTCH(FIN  ,KA  ,  KB  ,KCr  K  D  ,  KE  ,1  MP,  FINZ1  ,  .  .  . 

FINZ2 .FCUTZ1 ,FOUTZ2) 

PROCEDURAL 

IF(IMP.NE.     1.0)     GC   TO    10    ■ 

IF(KEEP.NE.  1.0)     GO    TO     10 

IF  (TIME. EC. 0.0)     GC    TO    10 

FOUT=KA*FIN+KB*FINZ1+KC*FINZ2-KD*FOOTZ1-KE*FOUTZ2 

FOUTZ2    =    FCUTZ1 

FOUTZ  1    =    FCUT 

FINZ2    =    FINZ1 

FINZ1     =    FIN 
10       CONTINUE 
ENDMAC 
* 

*  ***************************** 

*  *****FRECUENCY     AVERAGERS***** 

*  ***************************** 

MACRO      Z40=AV20U0  (Z20,IMP) 


PROCEDURAL 
IF  (KEEP 
IF     (TIME. NE. 0.0)     GO    TO    5 


IF     (KEEP. HE. 1 .0)    GO    TO    20 

(TIME. NE. 0.0) 
Z20Z1=Z20 


Z40Z  1  =  Z20 

DEL=0.0 

GO    TO     10 
5  IF     (IMP. EQ. 1.0)     GO    TC    10 

Z40=ZU0Z1+DEL 

GO    TC    15 
10  Z40=Z20 

DEL=  (Z2C-Z20Z1)  /2  .  0 

Z20Z  1  =  Z  2C 

Z40Z1=ZU0 
15  CONTINUE 

20  CONTINUE 
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ENDM 
MACR 


5 
10 


15 

20 

ENDM 
MACR 


10 
15 


AC 

C  Z80=A 
EROCSDUR 
IF  (KE 
IF  (TI 
Z40Z  1  = 
Z80Z1= 
DEL=0. 
GO  TC 
IF  (IM 
Z80  =  Z8 
GO  TO 
Z8  0  =  Z4 
DEL=  12 

Z40zi= 

Z80Z  1= 

CONTIN 
CO  NT  IN 

AC 

C    Z80=AV 

EBOCEDUR 
IF  (KE 
IF  (TI 
Z20Z  1  = 
Z80Z1= 
GO  TO 
IF  (IM 
Z80=Z8 
GO  TC 
Z80=Z2 
DEL2  =  ( 
CONTIN 
Z20Z1= 
Z80Z  1  = 
CONTIN 

AC 


V4080(Z40,IMP) 

AI 

EP.N 

BE.H 

zao 
zao 
o 

10 

P.EQ 
0Z1  + 
15 

0 

40-Z40Z1)  /2.0 

zao 

Z8C 

UE 
UE 


E.1.0)    GO    TO    20 
E.0.0)    GO    TO    5 


.  1.  0)     GO    TC    10 

DEL 


208C  (Z20,IMP) 

AI 

EP.N 

ME.N 

Z20 

Z20 

10 

E.EQ 

0Z1  + 

15 

0 

Z2C-Z20Z1 )/4. 0 

UE 

Z20 

Z80 

OE 


E.1.0)    GO    TO    20 
E.0.0)    GO    TO    5 


.1.0)     GO    TC    10 

DEL  2 


*************************************************** 


20 

ENDM 

* 

* 

*  *****INITIAL   CONDITIONS    AND    CONSTANTS    FOR    BA 

*  *****FLIGHT   CONDITION   OF    20,000    FT    AND    250    K 

*  *********** **  *  * ***************************** 

INITIAL 

CONSTANT  PS=97  2.4  9.QC=1 1 2 . 73, RAT  E20=0 . 05. R ATE4  0=0. 
PS3IC1=0. 0,PS3IC2=0. 0,PS4IC1=0.0 ,PR2IC1=0. 0,PR2IC2 
AA2IC2=0.0  ,NZ2IC1=0. 0,NZ2IC2=0.0 ,NZ 4Z 1=0. 0 .N Z5Z 1=0 
PK 10=0.46  4  7, FI6Z  1=0 . 0 .MP  1  1Z1  =  0. 0 , MP  1 1Z2=0  .  0 , MP1 2Z1 
P8A  =  0.6  90  84,E8E=-0. 9 9  06  8. E8C=0.  6 6312, P8D=- 0.99068, 
DS 1  =  0.0. DLE 1  =  0. 0,DTE 1=0. 0 ,LE2Z1  =  0. 0 , TE2Z1  =0.  O.LEDI 
LEDIC2  =  U.O  ,TEDIC1=0. 0 .TEDIC2=0.  0 , PR3Z 1=0.  O.PR4Z1=0 
RATE 8 0=0. 0  125, A A 2IC 1  =  0. 0,EK9=-1.  1543,MF12Z2=0.0,P8 

DYNAMIC 

SORT 
* 

*  **************** 

4c  *****!  jjpuTS***** 

*  **************** 

PS1=6.0*STEE(0.0) 

PR1=0.0 

AA1=22.0+RAMP  (4.0) 

NZ1=0.0 

RR1=0.0 

YR1=0.0 


*** 
*** 
*** 


,  •  .  . 

.0  ,. 
3539 


* 

* 

* 
* 
* 


************************************ 

*****FUNCTIONS     REQUIRING    NOSORT***** 
************************************ 
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*  ********************* 

*  *****FUNCTION     12***** 

*  ********************* 

RI=OC/PS 

F12T1=RI**2*9.6  25-.025*RI+  1.0 
F12T2  =  PS*7.96  9E-4+0  .  BH 
F12MAX=IIMIT{1 .0,8. 0.F12T2) 
F12T3=LIMIT  (1 .  0  ,F  1  2  MAX,  F1  2T  1) 
F1 2T5SLIHIT  (.  5    1.35    RI) 

F12T6=F12T5*(0!00  95  2*PS+4.  04)     +  (-PS*0. 0039  6-  1. 18) 
F12T4  =  LIMIT  (1  .  0  ,8.  0  ,F12T6) 
NCSORT 

IFj.NOT.RI.GT.0.5)     GO    TO    5 
F12    =    F12T4 
GO    TO    6 

5  CONTINUE 
F12    =    F12T3 

6  CONTINUE 
SORT 

* 

*  ********************* 

*  *****FUNCTION     40***** 

*  ********  4***** ***** ** 

F4  0T  1  =  3.  25-3.  0*  LI  MIT  (0.75,  0.8  5.  RI) 
F40T2=0.656  25-.0013  125*LIMIT(500.0.1300.0 .PS) 
F4  0T3= -0.2617  7+  (9.6  35E-4)  *  LIMIT  (500.  0,  180  0.0,  PS) 
PIQ=LIMIT(0.0, 1800. C. PS) 


F4  0T4=LIKIT  (0  '.  0  , 1  0  5  .  5,  QC-3  35.  0) 

F4  0T5=  (F40T4*  (-PIQ*  (1.T7428E-6)  +  (1.5  238 E-  2 

+0.475- (6 .5E-4) *PIQ 
F4  0T6=F4  0T3+(F4  0T2*LIMIT(0.7  5,0.85,RI) ) 
F40T7=F4CTU*(5.6746E-4+1.984lE-7*PIQ) 

+  J0.  16923-3.84615E-5*PIQ) 
F40T8=(C.  16  92  3-3.  86  4  15E-5*  PIQ)  ..  . 

+  (F4  0T4*(  1.  6  7  247E-3-9.2915  2E-7*PIQ)  ) 


NCSORT 

IF     (. NOT. RI.GE.  0.75)    GO    TO    20 
F40T9  =  F40r8*F40T1 
GC    TO    30 
20  CONTINUE 

F40T9=F40T8 
30  CONTINUE 

IF     (.NOT.  (QC.GE. 440. C. AND. RI.GE. 0.75))     GO    TO    40 
F40T10=F40T6 
GO    TO    50 
40  CONTINUE 

F40T  10  =  F40T7 
50  CONTINUE 

IF     (.NO!  .PIQ.GT.980.C)     GO    TO    60 
F40=F40T9 
GO    TO    80 
60  CONTINUE 

IF     (  .NOT.PIQ.LE  .500. C)     GO    TO    70 
F40=F40T5 
GC    TO    80 
70  CONTINUE 

F4OF40T10 
80  CONTINUE 

* 

*  ************************** 

*  *****SORTED   FUNCTIONS***** 

*  ************************** 

SCRT 

*  ************************* 

*  *****FUKCTICNS     22,23***** 

*  ************************* 

F22  =  0.0167*(-80  0.0  +  LIMIT(8  00.0,900.0,QC)  ) 
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F  23  =  3.  14  35-0.  14  29  *LIMIT  (15.0,21  .998,  ALPHAT) 

*  ******************************************* 

*  *****FUNCTIONS     25,27,28,29 U,29L,32A, 37***** 

*  **£************?*******************  ******** 

F25  =  47. 6  36-0.0  5106*LIMIT(6  00.0,835.  0,QC) 


* 
* 

* 

* 
* 


F27=  1.3  2  8*  (ALPHAT+ 7.  658  4-1 7.86*LIMIT (0.44 , 0. 63, R I) ) 
F2 8= 44. 5 5  1-0. 0  4 05 8* LIMIT (2 60. 0,9  50.  0 ,QC) 

F29U=87. 38 25-7 6. 25* LIMIT (0.7, 1. 146,  R I) 
F29L=0.0 


QKF=LIMIT  (200. 0,2  0  0  0.0, QC) 
F32A=1QC.0/QKF 

* 

F37  =  2.5-0. 5*LIMIT (3.C,5.0,  NZA) 

*  ********  ** *******  ********* 

*  *****FUNCTIONS    68,107***** 

*  ********  4* **********  ****** 

F6 8  =  -0.0  0  2977*  (-4  80.  OLIMI  T  (260  .  0,4  8  0. 0,QC)  ) 
* 

F  10  711=  (  (-5.71  4 E- 7)  *ABS  (QC -750.0)  )  +  (8.4E-4) 

F107  =  LIKIT  (0.  0, 10  00C.0,F10  7T1) 

* 

*  **************************  * 

*  *****IMPULSE   FUNCTIONS***** 

*  *************************** 

IMP20=IMPULS(0. 0,RA1E20) 
IMP40=IMEOLS(0. 0,RATE40) 
IMP8  0=IKPULS(O. 0,RATE80) 

* 

*  ********************** 

*  *****pK    FUNCTIONS***** 

*  ********************** 

PK1 1=1.65*F22 
PK12=0.5654 
PK16A=3.5*F32A 
PK17A=F37*F32A 
PK19A=F12*F32A*3.5 
PK21A=F12*F32A*F68 
PK22A=F12*F32A 
PK16=ZHCID  (IMP2  0,  PK16A) 
PK17=ZHCLD  (IMP20,PK17A) 
PK19  =  ZHCID  (IMP2  0,  PK19A) 
PK2  1=ZHOLD  (IMP20,PK21A) 
PK22  =  ZHCID  (IMP20,  PK22A) 
* 

*  ******************************** 

*  *****pjxCT    STICK    INIUT    PATH***** 

*  ******************************** 

PS 2= EE A DSP (-2. 0,2.0,PS1) 

PS3  =  CMPXPL  (PS3IC1  ,PS3IC2,0  .  14  ,  27.  3,  PS2*  106  .  4  7) 
PS4  =  REALPL  (PS4IC1, 7.S365E- 3,PS3) 
PS5  =  ZHCID  (IMP80,PS4) 
PS6=PS5 

PS7  =  PS6*  (7.0+  (0.2*AES(PS6)  )) 
PS8  =  LIMIT(-25.0,50.  0,PS7) 
PV3=FS8 

PV3A=ZHCLD  (IMP20,PV3) 
PS9=PV3*F32A 
* 

*  ****************************** 

*  *****piTCH    RATE   GYRC    PATH***** 
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* 
* 

* 


* 
* 


* 

* 
* 


****************************** 

PR2=CMPXFL  (PR2IC1,PE2IC2,0.89,7  8.5,PR1*6  16  2.  25) 
PR3=ZH0ID  (IMP80  ,PR2) 
P2A=1.0*PK  11*  (  1  .0-PK12) 


P2B=  H.0  +  PK11)  *  (1  .  0-EK12) 
P2C=1.0-FK12 


* 
* 
* 


PR4=ZLAG  (PR3.P2A,  P2  E ,P2C,  I  MP8  0  ,PR3Z  1  ,PR4Z  1) 
PV4  =  LIMIT(-80. 0.120.0,PR4) 
PR5=PV4*  <  (F68*F32A)  +  F40) 
PV4A=ZH0LD  (IMP20,  PV4) 

************************************* 

*****ANGIE    OF    ATTACK    SENSOR    PATH***** 
************************************* 

AA2=CMPXEL  ( AA 21 C1 , A A2IC2, 0 . 74 , 2  09. 0 , AA 1*4 368 1.0) 

AA3=ZHOID  (IMP4  0,AA2) 

ALPHAT=0.59*AA3+1.  9 

ALPHAS  =  AIPHAT 

AA6=ALPHAT-22. 0 

AA7=IIMIT  (0.0,  1000  0.0,  A  A6) 

PV1=PK17*AA7 

PV1  A  =  ZHCID  (IMP2  0,  PV1) 

*********************************** 

*****NOREAL    ACCELERCMETER    PATH***** 
*********************************** 

NZ2  =  CMPXFL  (NZ2IC1 , NZ2IC2, 0 . 89r 200.0 ,  NZ 1*4  000  0.0) 

NZ3  =  ZHOID  (IMP40,NZ2) 

NZA=NZ3 

NZ4=NZ3-  |RR1**2)*  (6.8529E-6) 

P5A=  (1.0  +  FK9*  (  1  .0-PK10)  ) 

P5B=  (1.0  +  FK9)  *  (1.0-EK10) 

P5C=1.0-PK  10 

NZ5  =  ZLAG  |NZ4.P5A,P5E,P5C,IMP4  0,NZ4Z1  ,NZ5Z1) 

PV2  =  LIMIT(-10.0,10.0,NZ5) 

NZ7=PV2*EK16 

PV2A=ZHOLD  (IM?20,PV2) 

********************************* 

*****FORWARD    INTEGRATOR    PATH***** 
********************************* 

FI2=MP3 

FT 3  =  11  MIT  (-10  000.  0-0.0,  FI 2) 

FI4=LIMIT(0.0, 10Q0  0.0,FI2) 

FI5=  (FI3*F23)  +  FI4 

FI6  =  ZINT(FI5,  0.  05  ,  1  .  C,I  MP2  0,FI6Z1) 

******************* 

*****jlATN    PATH***** 
******************* 

MP1  =  PV3A*PK22 

MP2=  (PV4A*FK21)  -MP1 

MP3=  (PK19*EV2A)  +PV1A+MP2 

MP4=LIMIT(-50.0  ,25.  C.FI6) 

MP5A=AV204  0  (MP4,IMP20) 

MP5  =  ZHOLD(IMP4  0,MP5A) 

MP6=NZ7+EF5 

HP7=MP6+FV1 

M?8A=MP7+  (  (RR1*YR1)  *F107) 

MP8=ZHOLD  (I MP 4  0, MP 8 A) 

MP9  =  AV4080  (MP8,IMP40) 

MP10=MP9+PR5 

MP1  1=LIMIT  (-25.0,25.C,MP10) 

MP12=ZNCTCH  (MP1  1,P8A,P8B,  P  8C,  P8D,P8E  ,  IMP80, .  .. 

MP11Z  1,MP1 1Z2,MP12Z1, MP12Z2) 
MP13=MF12-PS9 

************************* 
*****STAEILATOR    PATH***** 
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*  ************************* 

ST2  =  IIMIT(-24.0,10.  5, MP  13) 
ST3=ST2-DS 1 

ST4  =  LIMIT(-24.0,10.5,ST3) 
ST5  =  CNTZE  (0.0125.ST4) 

STADEF=CMPXPL (TEDIC  1  ,TEDIC 2 ,0 . 7, 40 .0 ,ST5* 16  0  0.0) 
* 

4c  ******************************** 

*  *****LEADING    EDGE    FLAP    PATH***** 

*  ******************************** 

LE2  =  ZINT  <ALPHAS,0.06  25,0.9  37  5,IMP80,LE2Z1) 

LE2T=ZHCLD  (IMP20,LE2) 

LE2T1A=IE2T*F27 

LE2T1  =  ZHCLD(IMP20.LE2T1A) 

LE  2  T  2=  LI  MIT  (0.0  ,F28.LE2T1) 

LE3  =  LIMIT  (C.0,F29U,LE2T2) 

*  ************************************** 

*  *****18    DEG/SEC    RATE    LIMIT,     20    HZ***** 

*  ************************************** 

NO  SORT 

IF  (TIME. NE. 0.0)     GO    TC    85 

LE4Z  1=LE3 
85  CONTINUE 

DEL  1  =  LEj-LE4Z1 

DEL  II *= LI  MIT (- 0.9, 0.9,  DEL  1 ) 

LE4A=LEUZ1+DELLIM 

IF     (IMP20.  NE.  1.  0)     GC   TO    90 

LE4Z  1=LEUA 
90  CONTINUE 

LE4  =  ZHOLD(IMP20  ,LE4A) 

*  *****  *************  ********************* 

SORT 

LE5  =  AV2080  (LE4,IMP20) 

DLE2P1  =  LE5  +  3.0 

DLE2=DEADSP  (-DLE2P1  ,ELE2P1  ,  DLE1) 

DLE3  =  DLE1  +  DLE2 

*  NOTE:POSITIVE    LEADING    EDGE    FLAP    PATH    ONLY 
LE6=LE5+DLE3 

LE7  =  IIMIT  (-3.  0,33.0,LE6) 

*  ************************************** 

*  *****18    DEG/SEC    RATE    LIMIT,     80    HZ***** 

*  ************************************** 

NOSORT 

IF (TIME.NE.O. 0)     GO    10    110 

LE8Z1=LE7 
110         CONTINUE 

DEL1=IE7-LE8Z1 

DEL L I £=LIM IT (-0.22 5,0.2  25, DEL1) 

LE8A=LE8Z1+DELLIM 

IF     (IMP80.NE.  1.  0)     GO    TO    120 

LE8Z1=LE8A 
120         CONTINUE 

LE8  =  ZHCID  (IMP8  0  ,LE8  J) 

*  *************************************** 

SORT 

LE9  =  LIMIl(-3.  0.  33.0,LE8) 
LE10=CNTZR  (0.  0125,LE9) 

LEFLAP=CMPXPL (LEDIC  1  ,LEDIC 2, 1 . 4 , 20. 0 ,LE 10*40  0.0) 
* 

*  ***************************** 

*  *****NOSCRT   FUNCTION    24    ***** 

*  ***************************** 

F24L1=22.53  8-20.51*LIMIT(0.27,0.66,RI) 
F24 L 2=32. 7 6-3 6. 0*LI MIT (0. 6 6, 0 . 9 1 ,RI) 
F24T1  =  LIMIT(0.  0,1  0000  .0  ,ALPHAT) 

F24T2=ALPHAT- ( 1 4. 8769-7 .69 23*LIMIT ( 0 .27, 0 . 9 1 ,RI)  ) 
F24 T 3=- 2. 0*LIMIT( 0.0,10  000 .0,F24T2) 
F24T4=1.4*  (F24T  1+F24T3) 
NOSORT 
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10 

15 

* 
* 
* 


IF     (. NOT. RI.GT.  0.66)    GO    TO     10 

F2UL  =  F24L2 
GO    TO     15 
CONTINUE 

F24L=F24L1 
CONTINUE 
F24  =  LIMIT(0.0,F24L,  F24T4) 


**** 
**** 
**** 

TE2  = 

TE2T 

TE2T 

TE3  = 

*  NOTE 

TE4  = 

TE5  = 

TE6  = 

TEFL 

TERMINAL 

METHOD    RKS 

TIMER    FINT 

PRINT 

LABEL    STAD 

LABEL    TIME 

PAGE    XYPLO 

OUTPU 

END 

STOP 

EN  D  J  C  E 


***************************** 

♦TRAILING    EDGE    FLAP    PATH***** 
***************************** 

ZINT  (ALPHAS  .0.03125,0.  096875 ,IMP 80 ,TE2Z1 ) 

=ZHCLD  (IMP20,TE2) 

1  =  TE2T*F24 

LIMIT  (0.0,F25,TE2T1) 

:POSITIVE    TRAILING    EDGE    FLAP    PATH   ONLY 

TE3+DTE1 

LIMIT  (-8.  0  .45.0.TE4) 

QNTZEfO.0125.TE5) 

AF=CKFXPL(TEDIC1,TEDIC2,1.4,20.0,TE6*40.0) 

FX 

IM=10.0,OUTDEL  =  0.0125, PRDEL=0. 1 250 , DELT=0. 01  25 

STAEEF,LEFLAP,TEFLAP 
EF    IN    DEGREES 

IN    SECONDS 

T 

T    TIKE,STADEF 
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APPENDIX    C 


AIR     CATA    SCHEDULES 

LABEI    FUNCTION    12 

INITIAL 

PARAMETER    PS=     (2 00 0 .  C,  1  00  0.  0  , 500.  0  ,  20 0  .  0) 

DYNAMIC 
S  C  RT 

RI  =  RAMP  (0.  0) 

F12T1  =  RI**2*9.6  25-.  025*RI+  1.0 

F12T2=PS*7.96  9E-4+0.84 

F12MAX=IIMITM  .  0,  8.  0.F1  2T2  ) 

F12T3=LIMIT  (1  .  0  ,F  1  2  MAX,  F1  2T1) 

F12T5  =  LIMIT(.5,  1.35.BI) 

F12T6=F12T5*fO. 00  95  2*PS+4. 04)     +  (-PS*0. 0  03  96- 1.18) 

F12T4=LIMIT  (1.  0,8.  0,F12T6) 
NCSORT 

0.5)     GO    TO    5 


IF(.NOT.RI.GT 
F12  =  F12T4 

GO  TO  6 

5 

CONTINUE 

F12  =  F12T3 

6 

CONTINUE 

* 

TERMINAL 

TIMER    FINTIM   =2.0,     CUTDEL=0.05,    PRDEL=0.05 
PRINT    F12T1,F12T2,F12T3,F1 2T4  ,F  1  2T5  ,  F1  2T6 ,F12MAX ,F12 
PAGE.XYPICT 
OUTPUT    RI,    F12 
END 
STOP 
SNDJCB 


LABEL    FUNCTION    24 

INITIAL 

PARAMETEF  RI=  (  0  .27 , 0.66 , 0.  91 ) 

* 

DYNAMIC 
* 

SORT 

ALPHA  =  RAMP  (0.  0) 

F24I 1=22.538-20 .51*LIMIT(0 .27,0.66, R I) 

F24L2=32.76-36. 0*LI MIT (0. 6 6 , 0 .9 1  ,RI) 

F24T1  =  LI  KIT  (0.0, 10000.0,  ALP  HA) 

F24T2=ALPHA-(14.87  6  9-7.692  3*LIMIT(0. 27,0. 91, RI)  ) 

F24T3=-2.0*LIMIT(0. 0,10  000 .0,F24T2) 

F24T4=1.4*  (F24T  1+F24T3) 

NOSORT 

IF    j. NOT. RI.GT. 0.66)     GO    TO     10 

F24L=F24L2 
GO    TO    15 
10  CONTINUE 

F24L  =  F24L1 
15  CONTINUE 

F2 4  =  11  MIT  (0.0,F24L,F24T4) 
* 

TERM INAL 

TIMER    FINTIM=3  0.0,OUTDEL=0 .  5.PRDEL  =  0 .5 

PRINT    F24L 1,F 2  4L2 , F 24T1 ,F2 4T2 ,F2  4T3 , F2  4T4 ,F2  4L, F 24 
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PAGE    XYEICT 
OUTPUT    ALPHA, F24 
END 
STOP 
ENDJCE 


LABEL    FUNCTION    29 

INITIAL 

PARAMETER    FS=  (1250.0) 

* 

DYNAMIC 

RI=RAMP  (0.0) 

F29U  =  87.38  25-7  6.25*LIMIT(0  .7,  1.  146,  HI) 
F29L=0.0 

« 

TERMINAL 

TIMER    FINTIM=1.5.     PRDEL=0.05,    OUTDEL=0.05 
PRINT    RI,F29U,F29L 
PAGE    XYFLOT 
OUTPUT    RI,F29U,F29L 
END 
STOP 
ENDJCE 


LABEL    FUNCTION    40 

* 

INITIAL 

PARAMETER    QC= ( 1 70  0. C,44 0. 0 ,33  5 .0) 
* 

DYNAMIC 
S  CRT 

PS=    0.0001    ♦   RAMP  (O.C) 

RI=QC/PS 

F40T  1=3. 25-3.  0*LIMI  T  (0.  75,  0  .  85  .RI) 

F4  0T2=0.65  625-.00  13  125*LIMIT  (500 . 0. 1 80  0 . 0  .PS) 

F40T3  =  -0.2  617  7+ J9.6  35E-4)  *  LI  MIT  (500.  0,  180  0.0,  PS) 

PIQ=LIMIT(0.0, 1800. C.PS) 


F4  0T4  =  LIKIT  (0  .  0  , 1  0  5  .  6,  QC-3  35.  0) 

"  i  (FaOT4*  (-PIQ*  (1.174T 

♦0.475-    6 .5E-4) *PIQ 


F4  0T5=  (FaOT4*  (-PIQ*  (1.T742  8E-6)  ♦  ( 1 .  5  238E-3)  )  ) 


F4  0T6=F4  0T3+(F4  0T2*LIMIT(0 . 75 ,0 . 85,  RI)  ) 
F4  0T7=F40T4*(5. 67  4  6E-4+ 1 . 9 841 E-7*PIQ) . .. 

♦  (0.  16923-3.  8  U6 15  E-5  *PIQ) 
F4  0T8=  (0.  16  92  3-3.  86  4  15E-5*  PIQ)  ... 

+  (F4  0T4*(  1.  67  247E-3-9.29152E-7*PIQ)  ) 
* 

NCSORT 

IF     (. NOT. RI.GE. 0.75)    GO    TO    20 
F40T9  =  F40T8*F40T1 
GO    TO    30 
20  CONTINUE 

F40TS=F40T8 
30  CONTINUE 

IF     (.NOT. (QC.GE. 440. C. AND. RI.GE. 0.75))     GO    TO    40 
F40T10=F40T6 
GO    TO    50 
40  CONTINUE 

F40T  10=F40T7 
50  CONTINUE 

IF     (. NOT. PIQ.GT  .980.0)     GO    TO    60 
F40=F40T9 
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GO  TO  80 
60     CONTINUE 

IF  7.NOT.PIQ.LE  .500. C)  GO  TO  70 

F40=F40T5 

GO  TO  80 
70     CONTINUE 

F4G  =  F40T10 
80     CONTINUE 

TERMINAL 

TIMER    FINTIM=2  000.0,PRDEL=20.  0,  OUTDEL=20.  0 
PRINT    PS,QC,RI,F4  0T7,F4  0T9  ,F4  0T10,F4  0 
PAGE    XYFIOT 
OUTPUT    PS,F40 
END 
STOP 
ENDJCE 


LABEL    GENERATED    FUNCTION 

*  THIS    FROGEAM  GENERATES    QC    FUNCTIONS    2  2,25,28,32a,     AND    68 

* 

INITIAL 

PARAMETEF    ES=1000.0 
* 

DYNAMIC 

QC  =  RAMP<0.0) 
* 

F2  2  =  0.0167*  (-8  0  0.0  + LIMIT  (8  00.  0,  900.0,QC)  ) 

F25= 4 7.6  36-0. 05 106*LIMI T (6 00. 0,8  35.0,QC) 

* 

F2  8  =  44. 551-0.  0  4  05  8*  LIMIT  (2  60.  0,9  50.  0,QC) 
* 

QKF=LIMIT (200. 0,2000.0,QC) 
F32A=100.0    /   QKF 

F6  8=-0. 0  02977* (-48  0 . 0+LIMI T  (260 . 0, 4  8 0. 0,QC) ) 
* 

TERMINAL 

TIMER    FINTIM=2500. 0.    OUTDEL=50.0.    PRDEL=50.0 

PRINT  QC  ,F22,F25,F28,QKF,F32A,F68 

PAGE  XYFLOT 

OUTFUT    QC,F22 

OUTPUT    QC,F25 

OUTFUT    CC,F28 

OUTPUT    QCF32A 

OUTFUT    CC,F68 
END 
STOP 
ENDJOE 


LABEI    GENERATED    FUNCTION 

INITIAL 

PARAMETER    RI=  (0 .44 , 0.63) 
* 
DYNAMIC 

ALPHA=    EAMF    (0.0) 

NZA=ALPHA 

ALPHAT=AIPHA 
* 

F23  =  3.  14  35-0.  14  29* LIMIT  (15  .0,  21  .998,  ALPHAT) 
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F27= 1.3  28*  (ALPHA* 7. 8  584-17 .86*LIMIT (0.44,0. 6  3, HI)) 
F37=2.5-0.5*L I MIT (3.0,5.0, NZA) 

TERMINAL 

TIMER    FINTIM=30.0,     PFDEL=1.0,    ODTDEL=1.0 

PRINT    F23,F27,F37 

PAGE    XYIIOT 

OUTPUT    ALPHAT,F23 

OUTPUT    AIPHA.F27 

OUTPUT    NZA,F37 
END 
STOP 
ENDJCE 
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APPENDIX    D 


COMPUTER    PROGRAMS    FOP.     SIGNAL    BLOCK    TESTING 


CEADSPACE    FUNCTION    AND   ALIASING    FILTERS 


INITIAL 

CONSTANT    PS2C1=2.0 
* 

DYNAMIC 

PS1=    -5.0    +    RAMP(O.O) 

PS2=DEADSP  (-FS2CT,PS2C1,PS1) 

PS3=CMPXPI  (0.0,0.0.0.  1U.2  7.3  ,PS2*106.4  7) 

FSU  =  REALPL (0.0,  7. 936  5E-3, PS3 ) 
* 


TERMINAL 

TIMER 

LABEL 

PRINT 

OUTPU 

PAGE    XYPLCT 

END 

SIOP 
ENDJCE 


FINTIM=10.0,    OUTEEL=0.25,    PRDEL=0.25 
STICK    DEADSPACE    PLOT 
P  c 1    PS  2 
TIME,    PS  1 ,PS2,PS3,PS4 
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FREQUENCY     AVERAGERS 

MACRO    Z80  =  AV2080  (Z20,IMP) 
PROCEDURAL 

IF     (TIME. NE. 0.0)     GO    TO    52 

Z20Z1=Z20 

Z80Z  1  =  Z20 

GO    TC    54 
52  IF     (IMP. EQ. 1.0)     GO    TO    54 

Z80=Z80Z1+DEL2 

GO    TO    56 
54  Z80  =  Z20 

DEL2=(Z2C-Z20Z1)/4.  C 
56  CONTINUE 

Z20Z  1  =  Z20 

Z80Z1=Z80 
ENDMAC 
INITIAL 

CONSTANT       RATE=0.05 
DYNAMIC 

Y1=    RAMP  (0.0) 

Y1 A  =  EXF  (Y1) 

Y2=    IMPUIS  (0.0, RATE) 

Y3=     ZHCID  (Y2,  Y1  A) 

Y4  =    AV2080(Y3,Y2 
TERMINAL 

METHOD  RKSFX 

TIMER    FINTIM=1.  0/OUTCEL=0.  0125,  PRDEL=0  .  0  1  25f  DELT  =0.  01  2  5 

LABEL    20    TC    80    AVERAGER 

PRINT    Y2,Y3,Y4 
END 
STOP 
ENDJOE 

MACRC      Z40=AV2040  (Z20  ,IMP) 
PROCEDURAL 

IF     (TIME. NE. 0.0)    GO    TO    5 

Z20Z  1  =  Z20 

Z40Z1=Z20 

DEL=0. 0 

GO    TC    10 
5  IF     (IMP.EQ.  1.  0)     GO    TO    10 

Z40=Z40Z1  +  DEL 

GO    TO     15 
10  Z40=Z20 

DEL=  {Z20-Z2OZ1)  /2.  0 

Z20Z1=Z20 

Z40Z  1  =  Z40 
15  CONTINUE 

ENDMAC 
INITIAL 

CONSTANT       RATE=0.05 
DYNAMIC 

Y1=    RAMP  (0.0) 

Y1A=EXP  (Y1) 

Y2=    IMPUIS  (0.  0,RATE) 

Y3=     ZHOLE  (Y2,  Y1  A) 

Y4  =  AV2040(Y3,  Y2) 
TERMINAL 

METHOD  RKSFX 

TIMER    FINTIM=1. 0,OUTDEL=0.  025 , PR DEL  =  0 . 025, DELT=0  .0125 

LABEL    20    TC    40     AVERAGER 

PRINT  Y1  ,YU,Y2,Y3,Y4 

PAGE  XYPLOT 

OUTPUT    TIME,Y1 A,Y3, Y4 
END 
STOP 
ENDJOE 
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BATE    LIMITERS 

INITIAL 
DYNAMIC 

Y1  =  RAMP  (0.0) 

Y2=EXP  HI) 

Y3  =  LIMIT  (CO, 3.  0,  Y2) 

YU=  (-1  .0+EXP( -TIME* 2.0)  )  *STEP  (2.0) 

Y5=Y3+Y4 

I MP 2 0=1 MPU IS (0.0,0. 05) 

IMPU0  =  IMPULS(  0.0,0.  025) 

IMP8  0=IMFULS(0.0, 0. 0125) 

Y20=ZHOLD(IMP20  ,Y5' 

Y40  =  ZHCID  (IMPUO  ,Y5 

Y80  =  ZHOID(IMP80,Y5! 
NOSORT 

IF  (TIME. NE. 0.0)     GO    TO    80 

LESZ1=Y20 
80  CONTINUE 

DEL1=Y20-IE8Z1 

DELLIK=LIMIT(-0.1 ,0.1,DEL1 ) 

LE8=IE8Z1+DELLIM 

IF     (IMP20.NE.  1  .  0)     GO    TO    90 

LE8Z1=LE8 
90  CONTINUE 

TERMINAL 

METHOD    RKSFX 

TIMER    FINTIM=5. 0,OUTDEL=0.  01 25, PR DEL =0 . 0 1 25, DSLT =0. 01 2 5 

PRINT    Y3,YU,Y5,DEL1  ,DELLIM ,LE8Z1  ,LE8  , Y20 , 1 MP 20 

PAGE    XYEICT 

OUTPUT    TIME,Y20,LE8 

OUTPUT    TIME,Y20,LE8Z1 
END 
SIOP 
ENDJCB 
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DIGITAL    FILTERS 

MACRO    FOUT=ZINT(FIN,  KA,KB  ,IMP,F0UTZ1) 

PROCEDURAL 

IF(IMP.NE.     1.0)      GO    TO     10 

IF  (KEFF.NE.  1.0)     GO    TO    10 

IF(TIME.EO.O.O)     GO    TO     10 

FCUT=KA*FIK    +    KB*FOUTZ1 

POUTZ  1=FOUT 
10       CONTINUE 
EN  DM  AC 
MACRO    FOUT=ZLAG  (FIN, KA, KB , KC, IMP , FINZ 1 , FOUTZ 1) 

PROCEDURAL 

IF(IMP.NE.     1.0)     GC    TO     10 

IF(KEEP.NE.  1.0)     GO    TO     10 

IF  (TIME.  EQ. 0.0)     GC    TO     10 

EOUT=KA*FIN    -    KB*FINZ1    ♦    KC*F0UTZ1 

FOUTZ1=FOUT 

FINZ1=FIN 
10       CONTINUE 
EN  DM  AC 

MACRO    FOUT=ZNCTCH  (FI  N  .K  A  ,  KE,KC,  K  E,  K  E  ,IMP  ,  .  .. 
FINZ1,FINZ2,FOUTZ1  ,  FOUTZ  2) 

EEOCEDURAI 

IE(IMP.NE.     1.0)     GO    TO     10 

IF  (KEEP.  NE.1  .0)     GC    TO    10 

IF(TIME.ECj.O.O)     GO    TO     10 

FOUT=KA*FIN+KB*FINZ1+KC*FINZ2-KD*FOJTZ 1-KE*FOUTZ2 

EOUTZ2    =    FOUTZ1 

FOUTZ1    =    FCUT 

EINZ2    =    FINZ1 

FINZ1     =    FIS 
10       CONTINUE 
ENDMAC 
INITIAL 

PA  RAM    INTN=.0125.INTD=1.0  ,AA11A1=0.0,AA11A=0.0 
PARAM    PK9=-1.15U3,    PK  10  =  0  .U647, A A1 0B1 =0. 0 , AA 1 1 B 1=0 .0 
PAR  AM    NA=. 6908  4. NB=-. 99 06 8,NC=. 6 6 31 2, ND=- . 99  06  8, NE  =  . 35 39  6 
PARAM    AA10C1=0.0,AA 1 0C2  =  0  .  C, A A1 1 C1 =0 . 0 ,AA 1 1C 2  =  0.0 
LA=  (1.+PK9*  (1  .-EK10)  ) 
LB=(  1.+PK9)*  (  1.-PK10) 
LC  =  T  .-PK10 
DYNAMIC 

AA9A    =    STEPjO.O) 
IMP=IMPULS  (0.  0,0.01  25) 
AA10A    =    ZHCLD  (IME,AA9  A) 

AA11A    =    ZINT  (AA10A,INTN,INTD,IMP  ,AA11A1) 
AA11  E=ZLAG  (AA10A,LA,LB,  LC  ,  IMP  ,  A  A  10B  1 ,  AA11B1) 
AA11C    =    ZNOTCH  (A  A  10  A , NA , N B  ,NC ,ND ,NE ,IMP, A A 10C1 ,... 

AA10C2,AA11C1,AA1  1C2) 
TERMINAL 
M P  "P  W  CD    RKSFY 

TIMER    FINTIM=1.5,OUTDEL=0.012  5,DELT=0.0125 
PRINT    AA10A,AA11 A,AA1  1B, AA11C 
OUTPUT    AA10A  (0.0,  1.  5)  ,AA  1  1A  (0  .  0,  1  .  5)  , .  . . 

AA11B  (0.0,1.5)  ,AA11C(0.0,  1.5) 
END 
STOP 
ENDJOE 
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